0^ 



p. 



X 



Theoretical Physics Institute 
University of Minnesota 



TPI-MINN-94/38-T 
;r ; UMN-TH-1320-94 

0\ ; November 1994 

> 

o 

Moments of lepton spectrum in B decays and the 

^ ■ rrib — rric quark mass difference. 

^' 

(N ■ M.B. Voloshin 

Theoretical Physics Institute, University of Minnesota 
Minneapohs, MN 55455 
(^ ■ and 

Institute of Theoretical and Experimental Physics 



§^; Moscow, 117259 



Abstract 

It is argued that the quark mass difference ruf, — rric can be extracted with a high accuracy 
from experimental data on ratia of moments of lepton energy spectrum in semileptonic 
decays of B mesons. Theoretical expressions for the moments are presented, which include 
perturbative as well as non-perturbative corrections. 



1 Introduction 

Masses of quarks as well as weak mixing angles are fundamental input parameters in the 
Standard Model. Thus it is important to know them with maximal possible accuracy. More- 
over, the values of the masses of the b and c quarks are correlated with the determination 
of the mixing parameter |Vbc| from the data on inclusive semileptonic B meson decay rates. 
Therefore an independent understanding of nib and rric is necessary for a precision determina- 
tion of I He I- There is a vast literature on extracting the values oirric and rrib from the data on 
charmonium, T resonances, charmed and B hadrons. However there is hardly a compelling 
argument in the literature to invalidate the original evaluations from the QCD sum rules of 
the 'on shell' quark masses: m, = 1.35 ± O.O5Gel^0Bi and m^ = 4.80 ± 0.03 GeV&, 
which can still serve as 'reference values' in discussion of dynamics of heavy hadrons. It 
is the purpose of this paper to point out that the accuracy of determining the difference 
between the quark masses m^ — rric can possibly be significantly improved by considering the 

ratia of the moments 

r dT 

M^ = jEr—dE, (1) 

of the lepton energy spectrum dT/dEi in semileptonic B decays. While theoretical expres- 
sions for the moments are sensitive to combinations of nib and nic (and also to \Vbc\) , their 
ratia for few first moments are to a high accuracy sensitive only to the mass difference 
nib — ""^c- Also it is natural to expect that experimentally the ratia of the moments are 
determined with better systematic accuracy, since the absolute normalization of the event 
rate cancels out in the ratia. 

On the theoretical side the ratia of the moments have the advantage of weak and 
controlled dependence on the infrared dynamics in QCD both perturbatively and non- 
perturbatively. For first few moments the perturbative corrections are expressed through 
as{nib) and the non-perturbative ones are suppressed by ?7i^^ and can be found by the OPE 
technique's B 111' H in terms of the quantities fx^/ml and f^l/ml with 

^^l = {B\{b7r%)\B) and ^,l = {B\ {b{cT ■ B)b) \B) , (2) 

where B is the chromo - magnetic field operator and tt = p — A is the covariant momentum 
operator for the heavy quark. The spin-dependent chromo - magnetic energy /i^ is related 
to the mass splitting of B* and B mesons: yU? = | {M^* — M%) k, 0.36 Gel^^, while for the 
kinetic energy fi^ only a lower bound exists^: ^^ > ^2 -^]-^ic]-^ follows from non-negativity 
of the operator (cr • tt)^ = tt^ — <j • B and an estimate'liy' /i^ ^ 0.5 ± 0.1 GeV'^ from QCD 
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sum rules. 

It should be also noticed that the difference mi, — rric, unlike each of the masses, is less 
sensitive to the infrared behavior in QCD and is a well defined quantity in QCD in the limit 
where both masses are heavy as compared to Aqqd- Indeed, because of confinement there 
is no real 'mass shell' for a quark. Therefore its mass can only be determined off shell and 
then extrapolated to a would-be on-shell value. For a heavy quark its mass can be found 
at a virtuality scale /i (i.e. at m'^ — p^ ~ 2 m fi) such that on one hand /i ^ Aqcd, which 
justifies a short- distance treatment, and on the other hand fi <^ m,. The latter condition 
ensures that the evolution of m(yu) towards the would-be mass shell does not depend on m 
in the leading order in 1/m. In particular, in the leading log approximation this evolution 
is described by the RG equationO' t3 

dm, 

where the constant c depends on the specific definition of the off-shell mass. The infrared 
singularity of a^ (infrared renormalon) prevents from integrating an equation like (§) down 
to /i = and thus really extrapolating the mass to the mass shell. However one can integrate 
the equation (|^) in any finite order in a^ and thus define the 'on-shell' mass of a heavy quark 
to a finite order. In this sense the 'on-shell' masses m^b and rric quoted above are the result of 
such extrapolation in the first order and are thus appropriate for using in other calculations 
in the first order in a^. Naturally this definition of quark mass changes with the order in 
as- However, increasing the order in a^ does not converge at a certain value of m because 
of the factorial divergence of the series in as caused by the infrared renormalon. A minimal 



residual error in the 'on shell' mass in this procedure is of the order of Aqcd ^ • On the 
other hand this uncertainty in a heavy quark mass does not depend on m in the limit of 
large m. Thus this uncertainty cancels in the difference of masses of two heavy quarks. As 
to the pre-asymptotic in the heavy quark mass limit corrections to the evolution equation 
(§), their contribution to the residual uncertainty is of the order of AqQjj/m,'^, which is quite 
small even for the charmed quark. 

The difference rrih — m^ can be estimated from the experimental values of the masses of 
D and B mesons: 

Ms-M:, = m,-m, + 4^ - 4^ + «K"'' "^^"') (4) 

Neglecting the terms, smaller than m~^ or m^^, and taking into account the inequality 



f^'i > f^l one finds a lower bound for the difference of tlie quark masses: 

mb-m^>MB-MD = 3.41 GeV . (5) 

Varying fi^ in the range from fi^ = fi"^ ^ 0.36 GeV^ up to fi^ = 0.6 GeV^ one finds nih — mc = 
3.44 ± 0.03 GeV, which is perfectly compatible with the quoted above estimates of each of 
the quark masses from QCD sum rules. 

Naturally, an independent measurement of this quark mass difference with a comparable 
or better accuracy would provide an additional consistency check for the heavy quark theory 
and, possibly, would enable a better quantitative understanding of the parameter /i^. As 
is discussed in the rest of this paper, a measurement of the ratia of few first moments (|I|) 
provides an excellent opportunity to independently determine the difference mb — rric. 

2 Moments of the lepton spectrum 

In the simplest approximation, where the QCD effects are neglected altogether the spectrum 
of charged lepton / in the decay b ^ clu is given by the well-known muon decay formula 

— = ToWo{x, n) (6) 

with To = Gl \Vcb\'^ ml/ {192 7v^) and 

wo{x, /i) = ^ "" ^}^ ~ ^ ~ ""^ [(1 - a:) (3 - 2x) + fi^3 - x)] , (7) 

(1 — xy L J 

where /i = mc/rrib and x = 2Ei/mb, so that the physical range of x goes from Xm = to 

xm = 1- /i^- 

With the first perturbative QCD correction and the first non-perturbative corrections, 
proportional to yU^/m^ and fP/m'^, taken into account the formula for the differential decay 
rate can be written as 

— -— = wo{x, /i) - 77— wi{x, n)^ 2 ^'rla;, At) + ^ Wg{x, /i) (8) 

L dx 6 IT mi m^ 

where the explicit expression for the perturbative correction function wi{x, fi) is extremely 
lengthy and can be found in the original papers [0, Q (see also [^, [1^)[[ The non- 



^ I am thankful to M. Jezabek, for pointing out to me the papers jlj] and ||lj|, where the calculation 
of the function wi{x, /i) has been finalized, and for sending me his and A. Czarnecki's FORTRAN code for 
numerical calculation of this function. 



perturbative correction functions Wt,{x, /i) and Wg{x, /x) are given bylif B 
3 (1 - xf 



w^{x, fi) = — — — 5 (-5-15/i^ + 20/i^ + 25a; + 21/i^a;-10/i6x-50x2- 



6 IJ,^ x^ + 2 ij,^ x^ + 50 x-^ - 25 X* + 5 x'') , (9) 

2x'^ (1 — u^ — x) 

Wg{x, /x) = ^ ^ — - {6 - 12 n^ - 30 fi^ - 13 X + 23 fi^ X + 20 ij,'^ X + 3 x^- 

3 (1 — a;) 

16 fj,'^ x^ - 5 fx^ x'^ + 9 x^ + 5 fx"^ x^ - 5 x^) . (10) 

The relative correction w-i (x, a)/wo(x, a) has a logarithmic singularity at the upper end- 
poi,. . = x„ of the spectrumB. wh.ch is a usual cousequeuce of the Sudakov form- 
factor liil. The relative non-perturbative corrections are still more singular: the ratio 
Wg{x, ^)/wo{x, /i) has a pole at x = xm and the ratio Wt^^x, fi)/wo{x, fi) has a double pole 
at the upper endpoint, which in particular reflects the difference in the kinematics of decay 
of a free heavy quark and of a heavy quark bound in hadronlii LLhI. This implies that the 
spectrum close to the endpoint is sensitive to the infrared hadron dynamics, while in integral 
quantities, like the moments of the lepton spectrum, the effects of this dynamics are inte- 
grated over and are present only in the form of small corrections. Naturally, this conclusion 
is valid only if the number n of the moment is not parametrically large, since high moments 
measure the spectrum near the upper endpoint, and all the infrared effects come back. In the 
expressions for the moments this growth of sensitivity to large distances reveals itself in the 
growth with n of the relative magnitude of the non-perturbative corrections. Therefore, one 
can consider as 'safe' the moments with such n, for which the non-perturbative correction 
is still small. As will be discussed, if one chooses to keep individual terms in the corrections 
at a level below 10% - 15%, this would limit the range of n to n < 5. Thus in what follows 
explicit results are presented for the moments in this range of n. 

According to eq.® ratio of the n-th moment to Mq (the total rate) r„ = Mn/M^ can be 
written as 

r„ = rf fl-?^e^ + 4e^ + 4^i^^l > (11) 



o - rt ■ 2 n 2 

3 TV mt mt 



where r^"-' is the same ratio in the lowest approximation: 

,(0) _ /"^fcV /o^"^ wo{x, iJ,) x'^ dx 



2 J J^ ^' wo{x, fi) dx 



and the corrections S!^' , S^^'' and S!;f' each being a function of /i are obtained from integrals 
with the corresponding correction function w{x, /i) in eq.(^ as 



'■1-m\„/'^ ,A^nj^ rl-M^ 



^ ^ /o ^ w{x, ij)x'^dx _ /o ^ w{x, ij)dx 

Jq~'^ Wo{x,fi)x^dx Jq~^ wq{x, ^)dx 

The moments of the lowest order function ^0(2;, /i) for n < 5 are listed in the Appendix. 

The correction coefficients 5^^^^ can in fact be found in a simple analytical form. This is 

possible due to the fact that the function w^(x, jj) is related to a modification by a small 

boost with < v^ >= fil/ml of the lepton spectrum described by the function Wo{x, /x) lif LL^: 

Integrating by parts one readily finds that 6^:^^ does not depend on n and is given by 

5W = n{n + 2)/6 . (15) 

The expression for the coefficients 5^^ can be found in a somewhat lengthy analytical 
form. The integrals in eq.(p!3D with the function Wg{x, fi) and n < 5 are listed in the 
Appendix. Similar integrals for the perturbative coefficients S!;^^ can also, perhaps, be done 
analytically as a function of the mass ratio /i. However, judging by the expression O' LLS 
for the function Wi{x, /i) and by the analytical expression for the dependence of the 0{as) 
correction to the total rate^c^, the resulting formulas should be prohibitively lengthy. For 
the practical purpose of analyzing experimental data it is sufficient however to have a table of 
these coefficients for values oi fi around the approximate actual value fi ~ 0.3. The numerical 
values of the coefficients S!;^^ and S^^ for n < 5 are given in Tables 1 and 2. Since for each 
n these coefficients are slowly varying functions of /i, the tabular values can be used for an 
interpolation. 

One can see from these numerical values and from eq.(0) that for a^ ~ 0.2, fil/ml ^ 
0.015 and fi'^/ml ^ 0.015 — 0.025 the perturbative correction to the ratia r„ is quite small 
as compared to the non-perturbative terms, and that each of the latter terms is within 10% 
- 15% range for n = 5, though the overall non-perturbative correction is significantly smaller 
due to a partial cancellation between the two terms. 



3 Discussion 

The estimates presented above illustrate that both the perturbative and the non-perturbative 
QCD corrections are sufficiently small and controllable in a number of ratia of moments of 



the lepton spectrum in semileptonic B decays, which number is sufficient for a detailed 
experimental study of the kinematical parameters of these decays. A simple numerical 
inspection reveals that the ratia r„, are in fact sensitive to the quark mass difference rrib — rric 
rather than to the individual quark masses. Therefore it is quite likely, that the value of 
this mass difference can be determined with high precision from experimental data, while 
to separate each of the masses, one will have to rely on other types of analyses, e.g. on the 
existing determination of rrib from the T sum rules, or, possibly, on one from the inclusive 
spectrum of photons, generated by the process 6 — > 57, which may become possible in a 



future development of the experiment ||2^ . 



One last remark is in order concerning a possible experimental measurement of the mo- 
ments Mn in e~^e~ annihilation at the T{4S) resonance. Since the resonance is slightly above 
the BB threshold, the B mesons have momentum of about 0.3 GeV, and their measured lep- 
ton spectrum is slightly distorted by boost. However in order to account for this boost in the 
integral quantities like the moments M„ there is no need to transform the measured lepton 
energy distribution to the B rest frame. The reason is that a small boost the expressions 
for the moments in the laboratory frame remain valid after adding in quadrature the 'intrin- 
sic' average momentum squared of the b quark in meson, /i^ with that of the B meson in 
the laboratory frame < p^ >. This obviously amounts to replacing the fi"^ by the effective 
quantity 

Jll = fil+<p'> . (16) 

One can notice that at the energy of the T(45') resonance the effect of the boost: < p^ >^ 
0.09GeV^, is rather small in comparison with /i^. 
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A Appendix 

The moments of the lowest-order energy distribution function, I^^ = Jq^^ wo{x, fj,)x^dx, 
entering eg. ([T^) , for n < 5 are given by the following expressions 

4°) = 1 - 8 /i^ + 8 ^'^ - / - 12 /i^ log(/i2) (17) 

/f ^ = ^ - ^ - 12/.^ + 20/.^ - ^ + ^ - 6;.^ (3 + /.^) log(/.^) (18) 

/f = ^ - ^ - 26/.^ + 32/.^ + i|^ - ^ - 8/ (3 + 2^^) log(/.^) (19) 

^"^°^ = ^-T/.^-^ + ^ + 5/ + /°-^ + /^-30/.^(l + /.^)log(/.^) (20^ 

5 48 u2 291 u^ 252 u^ „ ,. 12 u^^ 3ui^ , / , o\ , , 

14 7 5^5^^^^ ^^^35 75""^^^ (^ + ^^) ^°^^^ 

(0 



14 . „ „ .„ .„ 

(21) 

(0) _ 11 _ 27 /i^ _ 759 /x^ 329/^6 63 /x^ _ 7^ _ 7^ 9^ _ /i^ /x^ _ 
^ ~ 36 4 1(P^~6 ^~2 2 6~^^0 T^^" 

14 /i^ (3 + 5^2) log(^2) . (22) 

The first moments of the function Wg-. I^^^ = /o~^ Wg{x, fi) x^ dx, necessary for calcula- 
tion of the coefficients 6!^\ are given by the following expressions 

4^) = -^ + 4/.2-12/.^ + 12/.«-^-6/.Mog(/i2) (23) 

li'^ = -2 + ^ - 4/i^ + 8/i« - ^ + /i^° - 4/^2 log(;z2) (24) 

(25) 



6 6 

, 40 /i^ 28/i8 16/iio 5/xi2 / 10/ 

25 /i 4 u 2 — 3 u 

^^3 35 <d ^ \ ^ 3 

^(,) ^ _ 53 _ 4V 155^ _ ^ _ 25 u« + 8 «!" - ^^ + ^ 
^ 21 5 2 2 ^ '^ 30 14 

2 /i^ (^6 - 15 /i^ - 25 /i^) log(/i2) (26) 

(,) ^ _ 75 _ 76^ 1553^^ - 86 u*^ - ^^ + 20 u^' - ^-^^ + ^°^^'' - ^ - 
4 OS ^ in ^6 ^ 10 105 4 

2 /i^ (^8 - 27 /i^ - 60 /i^) log(/i2) (27) 



151 160 u^ 1299 u^ 1057 u^ 455 u« ,^ m 175 /i^ 
c;/i 7 c; Q 'i ^ Q 



15 14 27 



01 2 175 /^n 

'-21/i' ^j log(/i2) . 



103/.- 23;.-^5^_^^, , _2^^._^^P^g^^.^_ ^28) 
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rric/mb 


^!" 


4^^ 


4" 


4^^ 


^i" 


0.25 


0.0252 


0.0669 


0.1193 


0.1787 


0.2427 


0.26 


0.0226 


0.0619 


0.1120 


0.1693 


0.2312 


0.27 


0.0201 


0.0571 


0.1052 


0.1603 


0.2202 


0.28 


0.0178 


0.0527 


0.0987 


0.1519 


0.2098 


0.29 


0.0156 


0.0485 


0.0925 


0.1438 


0.1999 


0.30 


0.0135 


0.0446 


0.0867 


0.1362 


0.1904 


0.31 


0.0116 


0.0408 


0.0812 


0.1289 


0.1813 


0.32 


0.0099 


0.0373 


0.0760 


0.1219 


0.1726 


0.33 


0.0082 


0.0340 


0.0711 


0.1153 


0.1643 


0.34 


0.0066 


0.0309 


0.0663 


0.1090 


0.1563 


0.35 


0.0052 


0.0280 


0.0619 


0.1029 


0.1486 



Table 1: Numerical values of the perturbative correction coefficients 5^^ in eq.(ll) for n < 5 
and yU = mc/mi, in the range from 0.25 to 0.35 . 



rric/mb 


5[^^ 


5^^^ 


s!f^ 


5?^ 


4^) 


0.25 


-1.188 


-2.419 


-3.674 


-4.94 


-6.213 


0.26 


-1.187 


-2.414 


-3.664 


-4.926 


-6.195 


0.27 


-1.186 


-2.411 


-3.657 


-4.915 


-6.17 


0.28 


-1.185 


-2.408 


-3.651 


-4.905 


-6.165 


0.29 


-1.185 


-2.406 


-3.647 


-4.898 


-6.155 


0.30 


-1.185 


-2.405 


-3.644 


-4.893 


-6.147 


0.31 


-1.186 


-2.406 


-3.643 


-4.89 


-6.142 


0.32 


-1.187 


-2.407 


-3.644 


-4.89 


-6.14 


0.33 


-1.189 


-2.409 


-3.646 


-4.891 


-6.141 


0.34 


-1.192 


-2.413 


-3.65 


-4.895 


-6.144 


0.35 


-1.195 


-2.417 


-3.655 


-4.901 


-6.151 



Table 2: Numerical values of the non-perturbative correction coefficients 5^^^ in eq.(ll) for 
n < 5 and /x = rric/rnb in the range from 0.25 to 0.35 . 
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